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ABSTRACT 

Polarization of the > 400 keV hard tail of the microquasar Cygnus X-1 has been independently 
reported by INTEGRAL/IBIS, and INTEGRAL/SPl and interpreted as emission from a compact jet. 
These conclusions were, however, based on the accumulation of all INTEGRAL data regardless of the 
spectral state. We utilize additional INTEGRAL exposure accumulated until December 2012, and include 
the AMI/Ryle (15 GHz) radio data in our study. We separate the observations into hard, soft, and 
intermediate/transitional states and detect radio emission from a compact jet in hard and intermediate 
states, but not in the soft. The 10-400 keV INTEGRAL (JEM-X and IBIS) state resolved spectra are 
well modeled with thermal Comptonization and reflection components. We detect a hard tail in the 
0.4-2 MeV range for the hard state only. We extract the state dependent polarigrams of Cyg X-1, which 
all are compatible to no or undetectable level of polarization except in 400-2000 keV range in the hard 
state where the polarization fraction is 75 ± 32% and the polarization angle 40?0 ± 14?3. An upper 
limit on the 0.4-2 MeV soft state polarization fraction is 70%. Due to the short exposure, we obtain no 
meaningful constraint for the intermediate state. The likely detection of a > 400 keV polarized tail in 
the hard state, together with the simultaneous presence of a radio jet, reinforce the notion of a compact 
jet origin of the > 400 keV emission. 

Subject headings: accretion, accretion disks — black hole physics — stars: individual (Cyg X-1) — 
X-rays: stars 


1. INTRODUCTION 

Black hole binaries (BHB) transit through different 
“spectral states” during their outbursts. The two canon¬ 
ical ones are the soft state (also known as “high” state, 
hereafter HSS) and the hard state (also known as “low” 
state, hereafter LHS). 

In the HSS, the emission is dominated by a bright and 
warm (~ 1 keV) accretion disk, the level of variability is 
low, and the power density spectrum is power-law like. 
Little or no radio emission is detected in this state, which 
is interpreted as an evidence for an absence of jets. 

In the LHS, the disk is colder (< 0.5 keV) and might 
be trun cated at a large distance from the accretor (see 
however [Reynolds fc Millerl 120131 for examples of results 
contradicting this picture). The X-ray spectrum shows 
a strong power-law like component extending up to hun¬ 
dreds of keV, usually with a roll-over at typically 50- 
200 keV. The level of rapid variability is higher than in the 
HSS and the power density spectrum s shows a band lim¬ 


ited noise component, sometimes with quasi-periodic oscil¬ 
lations (QPOs) with frequencies in the range 0.1-10 Hz. 

A so- called “compact-jet” is detected mainl y through its 
emiss ion in the radio to infrared domain ('e.g.. lCorbel et ^ 
l2013f) where the spectrum can be modeled by oc 
where a < 0 up to a break fre quency Vhrpak that usu¬ 
ally lies in the inf rared domain ( Corbel &: Fendeil l2002t 
iRahoui et al.ll2011r) . Above this break, the spectrum is a 
typical synchrotron spectrum with a > 0.5. 

Intermediate states (IS) exist between the HSS 


and LH (see, e.g., 
iHoman &: Belloiiil l2005l : 


Remil lard fc McClintocld 120061 : 
Fender et al.l 120091 for 
and more detailed state classifications). Transitions be¬ 
tween states indicate drastic changes in the properties of 
the flows and thus are crucial for the understanding of 
accretion-ejection connection. During a transition from 
the LHS to the HSS, the conipact jet is quenched (e.g. 
iFender et al.lll999l:ICoriat et al.ll2Qll[) : discrete, sometimes 
superluminal ejections occur, the level of X-ray variability 
drops abruptly, and QPOs, if present in the given source, 
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change typ es before disappearing (see, e.g.. lVarniere et aP 
I 2 OIILI 2 OI 2 L and references therein for both a description of 
the different types of QPOs and a possible theoretical in¬ 
terpretation of their origin and link to the spectral states). 

The BHB Cygnus X-1 (Cyg X-1) was discovered in 
1964 (|Bowver et al.lflO^ . It is the first Galactic source 
known to host a black hole (|Boltonlll975ll as the primary, 
and recent esti mates led to a bla ck hole mass of Mbh = 
14.8 ± 1.0 Mq (lOrosz et all 1201111. The donor star is the 
09.7 lab star HDE 226868 (iBoltonlfl^ iWalbornlllOTl 
with Mhde = 19.2 ± 1.9 (lOrosz et all 1201 iP* . The 
systems orbital period is 5.6 d llWebster &: Murdin Il972t 
iBrocksopp et al.lll999l : iGies et al.ll2003H . it ha s an~Tnclina- 
tion i = 27? 1 ± 0?8 on the plane of the sky (|Orosz et al.l 
UMII), and i s located at a distance of d = 1.86 ± 0.12 kpc 
from Earth (|Reid et al.ll2011l : IXiang et al.ll2011ll . 

Cyg X-1 can be found in both the LHS and the HSS; up 
to 2010, it was predominantly in the LHS. Since then, the 
beha vior has changed and i t spends most of its time in the 
HSS (|Grinberg et al.l [201311 . It sometimes undergoes par¬ 
tial (“failed”) transitions from the LHS to the HSS, and 
can be found in a transit ional or intermediate state (e.g., 
iPottschmidt et aI1l2003bl). The detection of compact rel¬ 
ativistic jets in the LHS (|Stirling et aLll^OlHr) places Cyg 
X-1 in the family of microquasars. It is one of the few 
microquasars known to have a hard tail extending to (and 
beyond) the MeV range fe.g.. iMcConnell et al.ll2000h . 

The presence of the MeV tail has recently been con¬ 
firmed with the two main instruments onboard the 
ESA’s INTErnational Gamma-Ray Astrophysics Labo¬ 
ratory {INTEGRAL) -, the Imager o nBoard the INTE- 
GRALSatellite fIBIS: lUbertini et al.ll2003D and th e SPec - 
trometer on INTEGRAL ( SPI : IVe drenne et al. 2003fl . 
ICadolle Bel et al.l (|2006ll and lLa.urent et al. (l201lL : 
LRWll) detect it with IBIS, and iJourdain et M ( 20121 
; JRC12) with SPI. Utilizing photon that C ompton scat¬ 
tered in the upper plane of IBIS, (ISGRI; iLebrun et al.l 
[ 2 OOI. and a re ab sorbed in its lower plane PICsIT 
( Labanti et al.llToOSL sensitive in the 200keV-10MeV en¬ 
ergy band), we have shown that the > 400keV emission 
of Cyg X-1 was polarized at a level of about 70%. We 
obtained only a 20% upper limit on the degree of polar¬ 
ization at lower energies (LRWll). This result was in¬ 
dependently confirmed by JRC12 using SPI data. Both 
studies obtained compatible results for the properties of 
the polarized emission, reinforcing the genuineness of this 
discovery. Both teams also suggested that the polarized 
emission was due to synchrotron emission coming from a 
compact jet. 

The presence of polarized emission and the energy- 
dependency of the polarization have profound implication 
on our understanding of accretion and ejection processes. 
It can help to distinguish between the different proposed 
emitting media (Comptonization corona vs. synchrotron- 
self Compton jets) in microquasars and provide important 
clues to the composition, energetics and magnetic field of 
the jet. 

A problem of both studies, however, is that they accu¬ 
mulated all INTEGRAL data available to them, regardless 
of the spectral state and radio (and thus jet) properties of 
the source. In this paper, we separate the whole INTE- 


^®See also IZiolkowski II2014I 1 who obtain somewhat different values. 


GRAL data set (accumulated up to December 2012) into 
different spectral states and study the properties of the 
state resolved broad band 10-2000 keV emission of Cyg X- 
1. We also utilize the Ryle/AMI radio data to determine 
the state dependent level of 15 GHz jet emission. In §2, we 
describe the data reduction and in particular the so-called 
“Compton mode” that allowed us to discover and mea¬ 
sure the polarization in Cyg X-1 (LRWll). We carefully 
separate the data into spectal state following a procedure 
based based on the classification of iGrinberg et al.l (|2013ll 
that is described in (21 The results of both, the long term 
15 GHz radio monitoring and INTEGRAL state resolved 
analysis, are presented in §4. We discuss the implications 
of our analysis in §5 and summarize our results in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Standard data reduction of the JEM-X and 
IBIS/ISGRI data 

Cyg X-1 has been extensively observed by INTEGRAL 
since the launch of the satellite; prelim inary results of the 
very fir st observations are described bv IPottschmidt et al.l 
(I2n'03all . In LRWll, we first considered all uninter¬ 
rupted INTEGRAL pointings, also called “science win¬ 
dows” (ScWs), where Cyg X-1 had an offset angle of less 
than 10° from the center of the field of view. We then 
removed all ScWs from the performance verification phase 
so that our analysis covered the period from 2003 March 
24 (MJD 52722, satellite revolution [rev.] 54) to 2010 June 
26 (MJD 55369, rev. 938). We further excluded ScWs with 
less than 1000 s of good ISGRI and PICSiT time. This re¬ 
sulted in a total of 2098 ScWs. Polarization was detected 
when stacking this sample (LRWll). 

In the present work, we added to the LRWll sample 
all observations belonging to our Cyg X-1 INTEGRAL 
monitoring program (PI J. Wilms) made until 2012 De¬ 
cember 28 (MJD 56289.8, rev. 1246). We applied the 
same filtering criteria to select good ScWs. We consider 
both, the data of th e low energy detector of IBIS, ISGRI 
(Lebrun et al.| | 2003ll and of the X-ray monitors JEM-X 
( Lund et al.ll2003ll ~ISGRI is sensitive in the ^18-1000 keV 
energy range, although its response falls off rapidly above 
a few hundred keV. The two JEM-X units cover the soft X- 
ray (3-30 keV) band. All data were reduced with the Off 
Line Scientific Analysis (OSA) v.10. 0 software suite 
and the associated updated calibration jCaballero^^^ 
2 OI 2 II following the methods outlined by (iRodriguez et al. 

2008bD . 

The brightest and most active sources of the field, 
Cyg X-1, Cyg X-2, Cyg X-3, 3A 1954-b319, and EXO 
2032-1-375, were subsequently taken into account in the ex¬ 
traction of spectra and light curves. The Cyg X-1 spectra 
of each ScW were extracted with 67 spectral channels. All 
ScWs belonging to the same state were then averaged to 
produce one ISGRI spectrum for each spe ctral state. The 
spect ral classification is described in (ISlfsee lCrinberg et al.l 
I 2 OI 3 I for the details of the method). A systematic error of 
1% was applied for all spectral channels. Low significance 
channels were rebinned when necessary. 

The JEM-X telescopes have a much smaller field of view 
than IBIS. The number of ScWs during which Cyg X-1 can 
be seen by these instruments is therefore smaller. Since 
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rev. 983 (MJD 55501) both units are on for all observa¬ 
tions. Before rev. 983, both units were used alternately, 
so that the number of ScWs and time coverage is not the 
same. We derived images and spectra in the standard 
manner following the cookbook, separating the data into 
states following the same criteria as for the IBIS data. We 
applied systematic errors of 2% to all spectral channels of 
both units. 


2.2. The Compton Mode 


Thanks to its two position-sensitive detectors ISGRI 
and PICsIT, IBIS can be used as a Compton polarime- 
ter (|Lei et al.lll997ll . The concept behind a Compton po- 
larimeter utilizes the polarization dependence of the dif¬ 
ferential cross section for Compton scattering. 


da 


dfl 


E' 


En 


E' En 


^ = ^ —-b-^-2sin^6»cos^ 


( 1 ) 


^Eo E' 

where tq is the classical electron radius, E^ the energy of 
the incident photon, E' the energy of the scattered photon, 
6 the scattering angle, and (j) the azimuthal angle relative 
to the polarization direction. Linearly polarized photons 
scatter preferentially perpendicularly to the incident po¬ 
larization vector. Hence by examining the scattering angle 
distribution of the detected photons (also referred to as 
polarigrams) 

7V(()i) = S')!-I-aocos(2((() - ^o))] (2) 


where S is the mean count rate, one can derive the po¬ 
larization angle PA and polarization fraction H. With 
PA = (j)o — tt/ 2 + nn, equation [2] becomes = 

S[l —aocos(2(()) —PH))][27r]. The polarization angle there¬ 
fore corresponds to the minimum of N{<f)). The polar¬ 
ization fraction is H = ao/oioo, where oioo is the ampli¬ 
tude expected for a 100% polarized source, ob tained from 
Monte -Carlo simulations of the instrument. iForot et al.l 
( 20081 1 obtained oioo = 0.30 ± 0.02, hence a ~ 6.7% sys- 
tematics uncertainty that should be taken into account 
in the derivation of H. Recent simulations allowed us to 
reduce the uncertainty on oioo to ^ 3% which is small 
compared to the statistical errors on N{(f)) (see below), 
and is therefore neglected. 

To measure f V(d>). we followed the procedure described 
bv iForot et al.l (j2008ll that led to the successful detection 
of the polarized signal from the Crab nebula. We consider 
events that interacted once in ISCRI and once in PICsIT. 
These events are automatically selected on board through 
a time coincidence algorithm. The maximum allowed time 
window was set to 3.8 /iS during our observations. To 
derive the source flux as a function of (j), the Compton 
photons were accumulated in 6 angular bins, each with a 
width of 30° in azimuthal scattering angle. To improve the 
signal-to-noise ratio in each bin, we took advantage of the 
TT-symmetry of the differential cross section (Eq. [2|), i.e., 
the first bin contains the photons with 0° < cj) < 30° and 
180° < (j) < 210°, etc. Chance coincidences, i.e., photons 
interacting in both detectors, but not related to a Comp¬ 
ton event, were subtracted from each detector ima ge fol¬ 
lowing the procedure described bv iForot et al.l (l2008ll . The 
derived detector images were then deconvolved to obtain 
sky images. The flux of the source in each (()-bin was then 
measured by fitting the instrumental PSF to the source 
peak in the sky image. 


The uncertainty of is dominated by statistic fluc¬ 
tuations in our observations that are background domi¬ 
nated. Confidence intervals on oq and (po are not derived 
by a fit to the data, but thr ough a Bayesian ap¬ 

proach follo wing Forot et al.l (l2008ll . based on the work 
presented in IVaillancourfj ( 20061) . The applica bility of this 
metho d was recently thoroughly discussed in iMaier et al.l 
([20 i 3). In this computation, we suppose that all real po¬ 
larization angles and fractions have an uniform pr obabil¬ 
ity d i stribution (non-in formative prior densities; iQuinnI 
I 2 OI 2 I : iMaier et al.ll20l4) and that the real polarization an¬ 
gle and fraction are (po and uq. We then need the proba¬ 
bility density distribution of measuring a and p from Npt 
independent data points in N{(j)) over a period tt, based 
on Caus sian distributions for the orthogonal Stokes com¬ 
pone nts (IVaillancourtll2006l : iForot et al.l]200^ iMaier et al.l 

[20TI : 


dP{a, (p) = 


Npt 


exp 


NptS^ 

■ 2a| 


\a^ + Og — 2aao cos(2(;ii — 2(Pq)\ 


a da d(p 


(3) 


where as is the uncertainty of S. Credibility intervals of a 
and (p correspond to the intervals comprised between the 
minimum and maximum values of the parameter consid¬ 
ered in the two dimensional I a contour plot (not shown). 


3. THE SPECTRAL STATES OF CYG X-1 

We previously developed a method to classify the states 
of Cyg X-1 at any given time based on hardness and in¬ 
tensity measurements with the RXTE/ASM and MAXI 
and on the hard X-ra y flux obtained with the Swift/BAT 
(|Crinberg et al.l[2013l) . The ASM and MAXI data can be 
used to separate all three states (LHS, IS, HSS), while us¬ 
ing BAT we can only separate the HSS, but not the LHS 
and IS from each other. 

In microquasars in general and Cyg X-1 in particular, 
spec tral transitions ca n occur on short time scales (hours; 
e.g., iBock et n]l201lD . We thus performed the spectral 
classification for each individual ScW. This is a sufficiently 
small exposure over which to consider the source to be 
spectrally stable in the majority of the data. The IBIS 18- 
25 keV light curve is shown in Fig. [U where the different 
symbols and colors represent the three different spectral 
states. 

Where available, we used simultaneous ASM data for 
the state classification. ScWs with two simultaneous ASM 
measurements that are classified into different states are 
presumed to have ocurred during a state transition and 
excluded from the analysis. Most ScWs are, however, not 
strictly simultaneous with any ASM measurements. To 
classify ScWs without simultaneous ASM measurements, 
we therefore used the closest ASM measurement within 
6 h before or after a given ScW. For the remaining point¬ 
ings where no such ASM measurements exist, we used the 
same approach, first based on M AXI and then, i f neces - 
sary, based on BAT. As shown bv ICrinberg et al.l (|2013D . 
the probability that a state is wrongly assigned using this 
method is 5% within 6 h of an ASM pointing for both the 
LHS and HSS, and ^^25% for the IS. 
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Out of the 3302 IBIS ScWs, 1739 ('^52.7%) were taken 
during the LHS, 868 ('^26.3%) during the HSS and 316 
(^9.6%) during the IS. The remaining 379 (11.5%) have an 
uncertain classification and are thus not considered in our 
analysis: 351 lack a simultaneous or quasi-simultaneous 
all sky monitor measurements that would allow a classifi¬ 
cation and 28 were caught during a state transition. 

From MJD 55700 until MJD 56000, Cyg X-1 was 
highly variable and underwent several tran sitions on short 
times cales (particularly striking in Fig. 1 of iGrinberg et al.l 
1201411 . Inspection of some of the ScWs during this pe¬ 
riod shows that the source flux was very low in the ISGRl 
range. Since these ScWs were taken close to or in be¬ 
tween state transitions, they were also removed from our 
analysis in order to limit the potential uncertainties they 
could introduce into the final data products. This filtering 
removed 32 additional ScWs. 

The resulting IBIS exposure times are 2.05 Ms, 1.21 Ms, 
and 0.22 Ms for the LHS, HSS, and IS, respectively. The 
IBIS 18-25 keV and Ryle/AMI 15 GHz light curves are 
shown in Fig. [T] 

4. RESULTS 

4.1. Radio monitoring with the Ryle-AMI telescope 

The 2002-2012 15 GHz light curve classified utilizing 
the same method as for the INTEGRAL data is shown in 
Fig.m As known from previous studies of microquasars, 
including Gyg X-1, the LHS and the IS show a high level 
of radio activity, while the radio flux is at a level compat- 
ible with or very close to zero durin g the soft state (e.g., 
iGorbel et al.l 120031 : iGallo et 'n]l2003ll . 

Gyg X-1 shows a high level of radio variability with a 
mean flux (F'isghz) ~ 12 mJy and flares. The most promi¬ 
nent flare occurred on MJD 53055 (2004 February 20) 
and reached a flux of lldmJjO. Flares occur frequently 
(Fig.fl]) , and are connected with the X-ray beh avior of the 
source ([Fender et al.ll200^ IWilms et al.ll2007ll . The main 
flares seem to mostly coincide with the IS. The radio be¬ 
havior in the LHS is usually steadier (although there is, 
e.g., a flare at MJD 53700, Fig. [T]). 

We obtained state resolved mean radio fluxes of 
(Aisghz.lhs) = 13.5 mJy, (Fisghz.is) = 15.4 mJy, and 
{Ti 5 GHz,hss) = 4.6 mJy. Given the t ypical rms uncer¬ 
tainty of rms(5min) = 3.0 mJy (e.g. iPoolev &: Fended 
Il997t [Rodriguez et al.ll2008^ . the mean radio flux in the 
HSS is not detected at the 3a level. The few clear radio de¬ 
tections in this state (Fig.[T|) do not correspond to a steady 
level of e mission that would indicate a compact radio core 
(see also iFender et al.lIl999ll . Such detections preferably 
occur after radio flares and thus likely repre sent the relic 
emiss ion of previously ejected material (e.g. iGorbel et al.l 

I200I . 

4.2. State resolved spectral analysis 

We analyzed the state resolved e nergy spectra o f each 
instrument using XSPEC v. 12.8.0 (lArnaudl |l996f) . For 
JEM-X we considered the 10-20 keV rangel^ and 20- 


400 keV for the ISGRI data. We used Compton data be¬ 
tween 300 keV and 1 or 2 MeV, depending on their statis¬ 
tical quality. We started by fitting the spectra in the 10- 
400 keV range since the contribution of the 1 MeV power- 
law tail is supposed to be negligible here (LRWll). When 
an appropriate model was found for this restricted range, 
we added the 0.3-2 MeV Compton spectra and repeated 
the spectral analysis. 

In order to find the most appropriate models describ¬ 
ing the 10-400 keV spectra for all three spectral states, 
we proceeded in an incremental way. We started with a 
simple power-law model and added spectral model com¬ 
ponents until an acceptable fit according to y^-statistics 
was found. The phenomenological models were then re¬ 
placed by more physical Comptonization models. Since 
the > lOkeV data are not affected by the abs orption col¬ 
umn density of Nh y^a f ew 10^^ cm“^ ('e.g.. iBock et al.l 
l201ll : [&inberg et al.ll2015ll . no modeling of the foreground 
absorption is necessary. Cross-calibration uncertainties 
and source variability during different exposure times are 
taken into account by a normalization constant. The IS¬ 
GRl constant was frozen to one and the others left free 
to vary independently. We required that the values of the 
free constants were in the range [0.85,1.15]. The bestfit 
parameters obtained with the best phenomenological and 
with the Comptonization models are reported in Table |T] 
for all three spectral states. 

4.2.1. Hard state 

A power-law (Fig. [5^) gives a very poor representation 
of the 10-400 keV data with a reduced = xt ’^0 
for 87 degrees of freedom (dof). The broad band spec¬ 
trum clearly departs from a straight line, which indicates 
a curved spectrum. A power-law modified by a high energy 
cut-off (highecut*power in XSPEC) provides a highly sig¬ 
nificant improvement with = 0.89 for 85 dof (Fig. [2b)- 
We note, however, that the value of the high energy cut-off 
is unconstrained (Table |T]), and deviations are seen at high 
energies. 

As a cut-off power-law in this energy range is usually 
assumed to mimic the spectral shape produced by an in¬ 
verse Compton process, we replaced the p henomenological 
model by a more physical mod el comptt ( |Titarch uklll994 
iTitarchuk fc Lvubarskiil 119951 : iTitarchuk fc Huai Il995ll . 
Due to our comparably high lower energy bound, the tem¬ 
perature of the seed photons could not be constrained. W e 
fix it to a typical LHS value of 0.2 keV (I Wilms et al.ll2006ll . 
Although the high energy part is better represented, the fit 
is statistically unacceptable (y^ = 2.9 for 88 dof. Fig. |2j:) 
and significant residuals are visible in 10-80 keV range. In 
this range, hard X-ray photons can undergo reflection on 
the accretion disk. This effect is usually seen as an extra 
curvature, or bump in the 1 0-100 keV range in the spectra 
of microquasars and AGNs ([Garcia et al.ll20i3 . and refer¬ 
ences therein). We thus included a simple reflection model 
(reflect) convolved with the Gomptonization spectrum 
and fixed the inclination angle to 30°. The residuals are 
much smaller, yielding an acceptable fit (yj = 1.67 for 87 


^^This flare was studied in detail by IFender et all (120061 ). although the time axis of their Fig. 4 is wrongly associated with MJD 55049.2- 
55049.6 

^^Although the JEM-X units are well calibrated above ~5keV, we chose to ignore the lower energy bins. Due to the possible influence of a 
disk and presence of iron fluorescence line at 6.4 keV they would add uncertainties and degeneracies to the models that the low energy resolution 
does not allow us to constrain. 
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dof0. The normalization constants are 0.99, 1.06, and 
0.98 for JEM X-1, JEM X-2, and the Compton mode re¬ 
spectively. The best fit results of this model are listed in 
Table [U while the vFv broad band 10-400 keV spectrum 
is shown in the rightmost panel of Fig [21 

Next, we added the >400 keV Compton mode spec¬ 
trum (Fig. [3]) to our data. These points are significant 
at the 8.2 ct, 7.8ct, 6.9a, and 6.5tT levels in the 451.4- 
551.9keV, 551.9-706keV, 706-1000.8keV, and -l-2MeV 
energy bands. A significant deviation from the previous 
best fit model is clearly visible above 400 keV. An addi¬ 
tional power-law improves the fit significantly {xt = 1-4 
for 90 dof). The photon index of this additional power-law 
is hard (T = l.llg^) and marginally compatible with the 
values reported earlier (1.6 ± 0.2, LRWll). 

We model the reflection component using several dis¬ 
tinct approaches, assuming 1) that the power-law com¬ 
ponent undergoes the same reflection as the Comptoniza- 
tion component (ref lect* (power+comptt) ), 2) that both 
components undergo reflection but with different solid 
angles f2/27r (reflecti*(power)+reflect2*(comptt)), 
and 3) that the power-law is not reflected 
(power+ref lect*(comptt) ). In model 2, the value for the 
reflection of the power-law component is unconstrained. 
Model 3 results in a better xl, but the value of the power- 
law parameters are physically not acceptable. The nor¬ 
malization is too high and T too soft to reproduce the 
Compton data well. We therefore discard model 3, and 
consider model 1 as the most valid one. The value of the 
fit statistics is, however, still high with xt = 1-8 for 89 
dof. The residuals show deviations in the 50-200 keV 
range that may be due to the accumulation of a large 
set of ScWs taken at different calibration epochs and at 
different off-axis and roll angles of Cyg X-1. 

In the final fit to the 10keV-2MeV broad band spec¬ 
trum the parameters of the Comptonized component are 
kTf, = 53 ± 2keV and r = 1.15 ± 0.04, the reflection frac¬ 
tion is unchanged, and the photon index of the hard tail 
is r = 1.4 ^q3 , i.e., essentially compatible with the value 
reported in LRWll. The 0.4-1 MeV flux of the hard tail 
(Apow.0.4-1 MeV = 1-9 X 10“® erg cm“^ s“^) accounts for 
86 % of the total 0.4-1 MeV flux. 

4.2.2. Intermediate state 

A simple power-law gives a poor representation of the 
IS spectrum {xt ^ 28 for 82 dof). The inclusion of a high 
energy cut-off again greatly improves the fit (xt = 1-35, 
80 dof), but the normalization constant of the Compton 
spectrum (Ccompton) is outside the range that is consistent 
with the flux calibration of the instrument. Replacing the 
phenomenological model by a comptt model, with the seed 
photon temperature fixed at 0.8keVP^ leads to a good fit 
{xl = 1.1 for 81 dof) and a physical value of Ccompton- 
Although the statistics are not as good as in the LHS due 
to the much shorter total exposure, and the residuals in 
the 10-100 keV region are rather acceptable, we included a 


reflection component to be consistent with the LHS mod¬ 
eling. This approach improved the fit to xi = 0.97 for 
79 dof, which indicates a chance improvement of ^ 1.2% 
according to an F-test (the reflection component has a 
significance lower than 3 ct). The reflection fraction is very 
low and poorly constrained (H/27r = 0.04 ± 0.03). As the 
other parameters are not significantly affected, we still re¬ 
port the results obtained with the latter model in Table |T] 
despite this weak evidence for reflection. 

We then added the > 400 keV spectral points of the 
Compton mode spectrum (Fig. [3]). The previous model 
(without reflection) leads to a good representation {xt = 
1.2 for 83 dof), even if a deviation at high energies is 
present. Above 400 keV, however, only the 551.9-706keV 
Compton mode spectral point is significant at > Scr. 
Adding a power-law to the data leads to xl = 1-04 (81 
dof), which corresponds to a chance improvement of 0.1%. 
Not surprisingly, the power-law photon index L is poorly 
constrained (—0.4 < T < 1.8) and the normalization con¬ 
stant for the Compton mode spectral points Ccompton 
tends to a very low value if not forced to be above 0.85. 
Fixing r at 1.6 (LRWll, JRC12) and setting Ccompton 
to 0.90, the best value found in the LHS fits, allows us 
to estimate a 3o- upper limit for the 0.4-1 MeV flux of 
1.5 X 10“® ergcm“^ s“^. 

4.2.3. Soft State 

In the soft state the soft X-ray spectrum below lOkeV is 
dominated by a ^1 keV accretion disk. This disk, however, 
does not influence the >10keV spectrum studied here. A 
single power-law gives a poor representation of the broad 
band spectrum (y^ = 6.0 for 77 dof). The JEM-X and IS- 
GRI spectra are particularly discrepant and a possible hint 
for a high energy roll-over is seen in the residuals. A cut¬ 
off power-law greatly improves the fit statistics (y^ = 1.68 
for 75 dof), but the modeling is still not satisfactory and 
the JEM-X normalization constants are inconsistent with 
the detector calibration uncertainty. Residuals in the 15- 
50 keV range suggest that reflection on the accretion disk 
may occur. Adding the reflect model improves the fit 
to xt = 1-24 (74 dof). The normalization constant for 
the JEM-Xl detector has to be forced to be above 0.85 
as it would naturally converge to about 0.8. Alternatively 
describing the data with a simple reflected power-law with¬ 
out a cutoff results in xt = 1-88 (76 dof) and shows that 
the cut-off is genuinely present. The chance improvement 
from a reflected power-law spectrum to a reflected cut¬ 
off power-law spectrum is about 4 x 10“® according to 
an F-test. Replacing the phenomenological model with a 
Comptonization continuum that is modified by reflection 
yields xl = 1-69 (75 dof. Table |T]). 

The >400 keV data were then added. Since the Comp¬ 
ton mode spectrum has a very low statistical quality and in 
order not to increase the number of degeneracies, we fixed 
its constant to a value of one. The ref lect*comptt model 
is still statistically acceptable, and only the highest spec- 


have followed the recommendations of the IBIS user manual regarding the level of systematics added to the IBIS spectral points 
(http://isdc.unige.ch/integral/download/osa/doc/10.0/osa_umjbis/node74.html). Systematic errors of 1% seem to underestimate the real un¬ 
certainties when dealing with (large) data sets that encompass several calibration periods. We have tested this notion by refitting the spectra 
with 1.5% and 2% systematics and, as expected, obtain xi much closer to 1. Since the values of the spectral parameters, and thus the 
conclusions presented here, do not change significantly, we decided to present the spectral fits using the recommended 1% for systematics. 

^"^The choice of a higher seed photon temperature in the IS does not influence the > 10 keV spectrum and was only motivated by the desire 
to be consistent with the expectation of a hotter disk in the IS. 
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Table 1 

Spectral parameters of fit to the 10-400 keV spectra. Errors and limits on the spectral parameters are given 
AT THE 90% confidence LEVEL, WHILE THE ERRORS ON THE FLUXES ARE AT THE 68% LEVEL. 


reflect*highecut(powerlaw) 

State 

F 

Ecut 

(keV) 

Efold 

(keV) 

10-20keV 

Fluxes^ 
20-200 keV 

200-400keV 

LHS 

IS 

HSSf 

1.43 ±0.01 

1 Q7+0.02 

2.447 ±0.007 

< 12 
561^ 

155 ±4 
198 ±8 

4.48 ±0.02 

5.10 ±0.03 

1.83 ±0.01 

22.30 ±0.03 
18.47 ±0.04 

3.33 ±0.01 

3.56 ±0.03 
2.52 ±0.05 

0.18 ±0.03 

reflect*comptt 

State 

H/27r 

kTe 

T 


Fluxes^ 




(keV) 


10-20 keV 

20 - 200 keV 

200-400keV 

LHS 

0.13 ±0.02 

59.4^);° 

1.06 ±0.03 

4.37 ±0.03 

22.60 ±0.03 

3.70 ±0.03 

IS 

0.04 ±0.03 

54.4;3-6 

0.82 ±0.06 

5.30 ±0.03 

18.42 ±0.06 

1.93 ±0.07 

HSS 

0.36 ±0.03 

279 ± 15 

< 0.013 

1.84 ±0.09 

3.7 ±0.6 

0.36 ±0.01 


^ A reflection component with r2/27r=0.46±0.04 was included for a good spectral fit to be obtained 
1 In units of 10“® ergcm“^ s“^ 


tral bin (550-2000 keV) indicates an excess of source pho¬ 
tons compared to the model {xt = 1-76, 77 dof). Adding 
a powerlaw component to the overall spectrum improves 
the modeling of the spectrum to xt = 1-69 (75 dof), but 
if left free to vary F is completely unconstrained. We 
fix it at 1.6 {xt = 1'’70: 76 dof, 0.11 chance probability 
according to the F-test) and derive a 3 ct upper limit of 
0.93 X 10“® ergcm“^ s“^ for the 0.4-1 MeV flux. 

4.3. State resolved polarization analysis 

Polarigrams were obtained for each of the three spectral 
states. The low statistics of the IS (the source significance 
in the Compton mode is 5.70") does not permit us to con¬ 
strain the 250-3000 keV polarization fraction. 

In the HSS, the detection significances of Cyg X-1 with 
the Compton mode are 4.2cr and 7.7a in the 300-400 keV 
and 400-2000 keV bands respectively. A weighted least 
square fit procedure was used to fit the polarigrams. A 
constant represents the 400-2000 keV polarigram rather 
well {xt = 1-6, 5 dof, not shown). The la upper limit 
of the 400-2000 keV polarization fraction in the HSS is 
-70%. 

Figure |4] shows the LHS polarigrams in two energy 
ranges: 300-400 keV and 400-2000 keV. Note that the data 
selection described in Sec. |3] is at the origin of the differ¬ 
ent count rates between the polarigrams shown here and 
those in Fig. 2 or LRWll. This effect is more obvious 
at the highest energies where the different spectral slopes 
have the largest impact (Fig. |3l). In the lower energy 
range, the Compton mode detection significance of Cyg 
X-1 is 12.3(7. The polarigram can be described by a con¬ 
stant {xt = 0.87, 5 dof). We estimate an upper limit 
for the polarization fraction of 0300-400 keV.LHS = 22% in 
this energy range. The 400-2000 keV energy range polari¬ 
gram shows clear deviation from a constant and a constant 
poorly represents the data {xt = 3-5, 5 dof). The pattern 
and fit indicate that the photons scattered in ISGRI are 
not evenly distributed in angle, but that the scattering has 
a preferential angle as is expected from a polarized signal 


(§2.2). We describe the distribution of the polarigram with 
the expression of Eq.[5]and obtain the values of the param¬ 
eters oq and (j)o with a least square technique. The fit is 
better than that to a constant, although not formally good 
{xt = 2.0, 3 dof) due tothe low value of the periodogram 
at azimuthal angle 100° (Fig. 0]). Assuming that the non 
constancy of the polarigram is an evidence for a real po¬ 
larized signal (an assumption also consistent with the in¬ 
dependent detections_oFpolRrKed emission in Cygnus X-1 
with SPI; iJourdain et al.l[2M^ . we estimate that the po¬ 
larization fraction is > 10% at the 95% level, and greater 
than a few % at the 99.7% level. In other words we obtain 
a significance higher than 3a for a polarized 400-2000 keV 
LHS emission. We determine a polarization fraction and 
a credibility interval of n 4 oo_ 2 oookeV.LHS = 75 ± 32% at 
position angle PAcygx-i,LHS = 40.0° ± 14.3°. The val¬ 
ues of the polarization fraction and angle are consistent 
with those found by JRC12. We note that an error of 
180° in the angle convention in LRWll led to PA = 140°, 
which after correction becomes 40°, in agreement with our 
state resolved result (see also JRC12). Note that the good 
agreement of our results with those obtained with SPI in 
the case of Cyg X-1, but also the fact that IBIS also de¬ 
tected polarization in the high energy emission of the Crab 
nebula and pulsar (|Forot et al.l 20081) at values c ompatible 
with those obtained with SPI ( Dean et al.ll200^ . increase 
the confidence in the reality and validity of our results. 

5. DISCUSSION AND CONCLUSIONS 
5.1. The > 400 keV hard tail 

Our analysis shows that the > 400 keV spectral proper¬ 
ties of Cyg X-1 depend o n the spectral state, a s already 
suggesteil previously ('e.g.. lMcConnell et al.ll2003) . We de¬ 
tect a strong hard tail in the LHS, but obtain only upper 
limits in the IS and HSS. In Fig. [5] (left), we show a com¬ 
parison of our LHS 0.75-1 MeV spectrum with those ob¬ 
tained with CGRO/COMPTEL during the June 1996 HSS 
and d uring an extended LHS obse rved with this instru¬ 
ment (iMcConnell et al]l2002L [20001) . The IBIS hard tail is 
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compatible in terms of fluxes to within 2a to those seen 
with COMPTEL and also with the INTEGRAL/SVl hard 
tail reported in JRC12 (see their Fig. 5). It is worth not¬ 
ing that Cyg X-1 is a notoriously highly variable source 
and flux evolution between different epochs is expected. 
To illustrate this behavior, we study data taken during 
two main periods, the covering March 2003 to December 
2007 (MJD 52722-54460) and March 2008 to December 
2009 (MJD 54530-55196), respectively. These intervals 
correspond to the analysis of LRWll. The two Compton 
spectra are shown in Fig. [5] (right). It is obvious that the 
flux of the hard tail varied between the two epochs and an 
evolution of the slope may also be visible. We note that 
these data are not separated by states, but as they cover 
the period before December 2009, they are dominated by 
the LHS (e.g., Fig.[T|). We therefore conclude that, even in 
the same state, the har d tail shows luminosity variations. 

A major difference to lMcConnell et all (l2002 |l is the non - 
detection of a hard tail in the HSS. lMcConnell et al.l (|2002ll 
remark that between two epochs of HSS the hard tail ei¬ 
ther has a higher flux than in the ir averaged LHS or is not 
detected. iMalvshev et al.l (|2013f) also note that the 1996 
HSS hard tail comes from a single occurrence of this state 
and may not reflect a typical behavior. While detection 
of a hard tail during the LHS, even if it shows variations, 
indicates an underlying process that is rather stable over 
long period of times, the apparent transience of the HSS 
hard tail may indicate a different origin for this feature. 
For example, in Cyg X-3 the very high energy (Fermi) 
flares are associated with transitions to the ultra soft state 
()Fermi LAT Collaboration et al.ll2009tl . While this is sim¬ 
ilar to what is seen in Cyg X-1, all Fermi T'-ray flares of 
Cyg X-1 were reporte d during LHS (|Bodaghee et al.l[2013t 
IMalvshev et al.l[2013ll . While the lack of simultaneous very 
high energy GeV-TeV data, and/or polarimetric studies of 
the Comptel HSS hard tail prevent any further conclusion 
on the origin of this component, we conclude that the HSS 
hard tail and the LHS hard tail could be of different origin. 

5.2. Polarization of the hard tail and its possible origin 

A clear polarized signal is detected from Cyg X-1 with 
both IBIS and SPI while accumulating all data regard¬ 
less of the spectral state of the source (LRWll, JRC12). 
Here we separated the data into different spectral states 
and studied the state dependent high energy polarization 
properties. The HSS and 300-400 keV polarigram can be 
described by a constant, indicating no or very little levels 
of polarization. The 400-2000 keV LHS polarigram shows 
a large deviation from a constant, which indicates that the 
Compton scattering from ISGRI is not an evenly circular 
distribution on PICSiT. We interpret this as evidence for 
the presence of polarization in this energy range, even if 
the polarigram shows some deviations from the theoreti¬ 
cally expected curve (Fig. 0]). In the 300-400keV range, 
no evidence for polarization is found. This is consistent 
with the results obtained with SPI (JRC12) and shows 
that the polarization fraction is strongly energy depen¬ 
dent. Assuming a non-null level of polarization, a reason¬ 
able assumption given the above arguments, we estimate 
that the detection of polarized emission is significant at 
higher than Za. Given the larger amount of data and the 
separation into different spectral states, one could have 


expected to obtain more robust results, and more con¬ 
strained polarization parameters compared to LRWll and 
JGR12, while the uncertainties we obtained are, at best, 
of the same order. The reason for this lies in the behav¬ 
ior of Gyg X-1. While the earlier studies did not perform 
a state dependent analysis, our state classification reveals 
that both studies were dominated by data taken during 
the LHS and the IS (see also Fig. [T]). In fact, compared 
to LRWll the number of ScWs measured during the LHS 
increased only by 13%, and 75% of the IS data comes from 
the period considered in these earlier studies, i.e., before 
^MJD 55200. On the other hand, about 94% of the HSS 
data are new. It is therefore not surprising that our hard 
state result is consistent with the earlier results, as these 
were fully dominated by the (polarized) hard state data. 
The exposure of the less polarized IS data was not large 
enough to significantly “dilute” the polarization. Our re¬ 
sult clearly shows that larger amounts of data in each of 
the states are needed to refine the constrain on the polar¬ 
ization properties and their relation to the spectral states. 

Our LHS spectral analysis shows that the 10-1000 keV 
spectrum can be decomposed into several spectral com¬ 
ponents. A (reflected) thermal Gomptonization compo¬ 
nent in the range ^ 10-400 keV, and a power-law tail 
dominating above ^400keV f M.2.11) . The origin of the 
seed photons for the Gompton component may either be 
an accretion disk or synchrotron photons from a compact 
jet undergoing synchrotro n self-Comptonizatio n (SSC). As 
discussed elsewhere fe.g.. iLaurent et aLll2011tl . due to the 
large number of scatterings the Gomptonization spectrum 
from a medium with an optical depth t > 1 is not ex¬ 
pected to show intrinsic polarization, especially with such 
a high fraction as the one we detect here. Multiple Comp¬ 
ton scattering, as expected in such a medium, will “wash 
out” the polarization of the incident photons even if the 
seed photons are polarized (e.g., jet photons). It is there¬ 
fore not surprising that the LHS <400 keV component 
has a very low (or zero ) polarization fraction (see also 
iRussell fc Shahb^l2013H . 

The origin of the hard-MeV tail is much less clear and 
is highly model-dependent. Two main families of models 
can explain the (co-)existence of cut-off power-law and/or 
pure power-law like emissions in the energy spectra of 
BHBs. The first is based on the presence of an hybrid 
thermal-non thermal population of electron in a “corona”, 
and hybrid Gomptonization have been successfully applied 
to the spectra of, e.g., GRO J1655—40 or GX 3 39—4 (e.g., 
ICaballero-Garcfa et al.l[200^ iJoinet et al.ll2007ll , and even 
to ^ 1-1 0000 keV spectra of Gy g X-1, in both the HSS 
and LHS (|McGonneli et al1l2002li . A double thermal com¬ 
ponent has been used to r epresent the 20-1000 keV SPI 
data of IE 1740—2947 well (iBouchet et al.l[200^ . 

The second family of models essentially contains the 
same radiation processes but is based on the presence of a 
jet emitting through direct synchrotron emission, thermal 
Gomptonization of the soft disk photons on the base of the 
jet, and/or SS C of the synchrotron photons by the jet’s 
electrons (e.g., iMarkoff et al.l [200^ . Such models have 
been used with success to fit the data of som e BH, and par¬ 
ticula rly GX 339— 4 and XTE J1118-I-480 (|Markoff et al.l 
120051 : iMaitra et aDl2009l: iNowak et'n]l201lH . It should be 
noted, however, that they obviously rely on the presence 
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of a compact jet and therefore can only be valid when/if 
a compact jet is present. 

The capacities of the current high energy missions, al¬ 
beit excellent, do not permit today an easy discrimination 
between the different families of models. This has nicely 
been illustrated in the specific case of Cyg X-1 through 
the 0.8-300 k eV spectral analys is of the Suzaku/RXTE 
observations (|Nowak et al.l 1201111 . Different and comple¬ 
mentary approaches and arguments can, however, be used 
to try and identify the most likely origin for the hard 
power-law tail. The fundamental pla ne between radio lu- 
minosity and X-ray luminosity (e.g., ICorbel et al.l 1200^ 
has originally been used as an argument for a (small) con¬ 
tribution of the jet to the X-ray domain. Our approach 
here was to first separate the radio data sets into different 
spectral states, based on a model-independent classifica¬ 
tion. It was only because of this separation that could 
we show that a hard polarized tail is present in the LHS, 
with a high polarization fraction, together with a radio 
jet. The very high level of polarization can only orig¬ 
inate from optically thin synchrotron emission (see be¬ 
low), and implies a highly ordered magnetic field similar 
to that expected in jets. The synchrotron tail coincides 
with thermal Comptonization at lower energies and with 
a radio jet. The presence of all these features is qualita¬ 
tively compatible with the multi-component emission pro- 
cesses predict e d by compact jet models such as those of 
iMarkoff et al.l (120051) . While polarization of the (compact) 
jet emiss ion has already been w idely reported in the radio 
domain (|BrocksoDD et al.ll2013L and references therein), a 
thorough study of the multi-wavelength polarization prop- 
erties of Cyg X-1 h as on ly recently been undertaken by 
iRussell fc Shahb^ l|2013[) . These authors, in particular, 
showed that the multi-wavelength spectrum and polariza¬ 
tion properties of Cyg X-1 are quantitatively compatible 
with the presence of a compact jet that dominates the 
radio to IR domain and is also responsi ble for the MeV 
tail. A similar conclusion was drawn by iMalvshev et al.l 
(|2013f) under the assumption that “the polarization mea¬ 
surements were robust” and that therefore the MeV tail 
could only originate from synchrotron emission. This pre¬ 
requisite is reinforced by our refined and spectral-state de¬ 
pendent study of the polarized properties of Cyg X-1 at 
0.3-1 MeV. A polarized signal is indeed suggested by the 
400-2000 keV LHS data and the large degree of polariza¬ 
tion measured in ( jd.dl imolies a very ordered magnetic field. 
Assuming the magnetic field lines are anchored in the disk 
implies that the 7-ray polarized emission comes from close 
to the inner ridge of the disk where the magnetic energy 
density is highest. This constraint is difficult to reconcile 
with a spherically shaped corona medium, where the mag¬ 
netic field lines are more likely to be tangled and where 
the polarized component necessarily comes from the outer 
shells of the corona where the optical depth is the lowest. 
We, therefore, conclude that in the LHS the 0.4-2 MeV 
tail detected with INTEGRAL/I^IS is very likely due to 
optically thin synchrotron emission, and that this emission 
comes from the detected compact jet. 


5.3. Compatibility of the Cyg X-1 hard state parameters 
with synchrotron emission 


The synchrotron spectrum of a population of electrons 
following a power-law distribution, dN{E) oc E~Pd{E), 
where p is the particle distribution index, can be ap¬ 
proximated by a power-law function over a limited 
range of energy, i.e., Fp] p]^(E) oc E~°' over [Ei,E 2 ] 
(|Rvbicki &: Lightmanlll979l) . In the case of compact jets, 
two main domains are usually considered: the optically 
thick regime with a < 0 and the optically thin regime 
with a > 0. The break between the two regimes is typ¬ 
ically in the infrared band (e.g., ICprbelJ^Jjendea 2002|; 
Rahoui et al.ll20lit IRussell fc Shahbazjl2013HCorbel et al l 


2OI3I) . The iet synchrotron spectrum is optically thin to¬ 
wards higher energies, and here a = (j> — l)/2. The jet 
emission becomes negligible in the range from the opti¬ 
cal to hard X-rays when compared to the contribution 
from other components such as the companion star, the 
accretion disk, or the corona. Our results indicate that 
it again dominates abo ve a few hundred keV (see also 
IRussell h Shahba'3120131 for a multi-components model¬ 
ing of the multi-wavelength Cyg X-1 spectrum). 

The degree of polarized emission expected in the op- 
tically thin regime is H = (|Rvbicki &: LightmanI 

119791) . Since the flux spectral index a = L — 1, where V is 
the photon spectral index, it is possible to deduce p from 
the measured spectral shape. With L = l.d^oJ f H4.2.1l) we 
obtain p = l.Stg;^ and thus expect Hexpected = 

This value is consistent with the value measured in the 


LHS (H = 75 ±32%). 

We determine a polarization angle of ^40° compatible 
with the SPI results (JRC12). This angle, however, dif¬ 
fers by about 60° from the polarization angles m easured 
in the radio and optical (IRussell fc Shahb^l2013L and ref¬ 
erences therein). Assuming that the 7-ray polarized com¬ 
ponent comes from the region of the jet that is closest to 
the launch site (i.e., the inner ridge of the disk and/or the 
BH), the change in polarization angle implies that the field 
lines in the 7-ray emitting region have a different orienta¬ 
tion than further out in the jet. This could be indicative 
of a relatively large ope ning angle at the basis o f the jet 
(as sketched in Fig. 3 of IRussell fc Shahbadl2013l) . Alter¬ 
natively the large offset could simply be due to strongly 
twisted magnetic field lines close to the accretion disk. 


6. SUMMARY 

We have presented a broad band 10keV-2MeV spec¬ 
tral analysis of the microquasar Cyg X-1 based on about 
10 years of data collected with the JEM-X and IBIS tele¬ 
scopes onboard the INTEGRA L observatory. We h ave 
used the classification criteria of iGrinberg et al.l (120131) to 
separate the data into hard (LHS), intermediate (iS), and 
soft (HSS) states. We have studied the radio behavior of 
the source associated with these states as observed with 
the Ryle/AMI radio telescope at 15 GHz. 


• The <400 keV emission is well represented by 
Gomptonization spectra with further reflection 
off the accretion disk. In the LHS and IS, the 
Gomptonization process is thermal and optically 
thick, i.e., the corona has an optical depth r > 1 
and the electrons have a Maxwellian distribution. 
As a result, a clear high-energy cut-off is seen in 
the spectra. In the HSS the situation is not that 
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clear, although a cut-off is also necessary to give a 
good representation of the broad band spectrum. 

• A clear hard tail is detected in the LHS when 
also considering the 0.4-2 MeV data. This high 
energy component is well represented by a hard 
power-law with no obvious cut-off. The detection 
of the hard tail is compatible with earlier claims 
of the presence of such a component in spectra of 
Cyg X-1. We show that this component is variable 
within same state, as is seen when considering 
INTEGRAL observations from two arbitrary 
epochs but the same state. 


• In the radio domain, the 15 GHz data show a 
definite detection with averaged flux densities of 
respectively ^13 and ~15mJy in the LHS and 
IS, compatible with the presence of a compact jet 
in those states. No persistent radio emission is 
detected in the HSS, implying the absence of a 
compact radio core. 


• In the LHS, we measure a polarized signal 
above 400 keV with a large polarization fraction 
(75 ± 32%). This high degree of polarization and 
the polarization angle (40?0 ± 14?3) are both 
compatible with previous studies by us and others. 
We obtain non-constraining upper limits on the 
polarization fraction in the IS and HSS, which 
have significantly lower exposure. 


• The high degree of polarization of the hard tail 
can only originate from synchrotron emission in an 
highly ordered magnetic field. The demonstrated 
presence of radio emission in the LHS points 
towards the compact jet as the origin for the 0.4- 
2 MeV emission, corroborating earlier theoretica l 
and multi-wavelengths studies (jMalvshev et al.l 
l2013t [Russell &: Shahbad [201311 . Our spectral 
state-resolved and multi-wavelength approach 
therefore further confirms the conclusion presented 
in the earlier, non state-resolved studies based 
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Fig. 1.— The two plots represent the long term radio (15 GHz, upper panels) INTEGRAL/IQIS hard X-ray (18-25 keV) light curves (lower 
panels) of C yg X-1. The different s ymbols and colors (online version) indicate the different states of the source according to the classification 
presented bv IGrinberg et ^ <1201311 . The upper plot covers the period MJD 52760—54250, when the radio monitoring was made with the older 
implementation of the radio telescope (known as the Ryle telescope). The lower one covers MJD 54500-56800 and shows the radio monitoring 
made after the upgrade (AMI telescope). Note the different vertical scales for the radio fluxes between the two plots. 
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Fig. 2.— Left: residuals to the hard state 10-400 keV spectra. From a) to d) simple power-law, power-law with exponential cut-off, 

Comptonization, Comptonization convolved by a reflection model. Right: 10-500 keV v — Fu (JEM-X-I-ISGRI) spectrum with the best 
model (Comptonization convolved by a reflection model) superimposed as a continuous line. For the sake of clarity only one JEM X spectrum 
is represented. 
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Fig. 3.— Unfolded lOkeV-2 MeV INTEGRAL spectra of Cyg X-1 in the three spectral states. The LHS spectrum is in blue circle, the IS 
in green diamonds and the HSS one in red triangles (colored points are available in the online version). 
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Fig. 4.— Polarigrams obtained in the LHS. Left: 300-400 keV. Right: 400-2000 keV. 
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Fig. 5. — Left: Comparison of the INTEGRAL/YQIS and CGi20/Comptel high energy spectra. Right: Comparison of the INTE¬ 
GRAL/Qompton spectra accumulated over two different epochs. 
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